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Abstract. Motivated by a recent generalization of the Balian-Low theorem and by new research 
in wireless communications we analyze the construction of Wilson bases for general time-frequency 
lattices. We show that orthonormal Wilson bases for L 2 (R) can be constructed for any time- frequency 
lattice whose volume is |. We then focus on the spaces and C L which are the preferred settings 

for numerical and practical purposes. We demonstrate that with a properly adapted definition of 
Wilson bases the construction of orthonormal Wilson bases for general time-frequency lattices also 
holds true in these discrete settings. In our analysis we make use of certain metaplectic transforms. 
Finally we discuss some practical consequences of our theoretical findings. 
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1. Introduction. Gabor systems have become a popular tool, both in theory 
and in applications, e.g., see EOl El- However one drawback is that due to 
the Balian-Low theorem it is impossible to construct (orthogonal) Gabor bases for 
L 2 (R) with good time- frequency localization [201 • In |35U11| it has been shown that a 
modification of Gabor bases, so-called Wilson bases provide a means to circumvent the 
Balian-Low theorem. Indeed, there exist orthogonal Wilson bases for L 2 (R) whose 
basis functions have exponential decay in time and frequency. These Wilson bases 
can be constructed from certain tight Gabor frames with redundancy 2. 

Gabor frames are usually associated with rectangular time-frequency lattices, but 
they can also be defined for general non-separable lattices, see e.g. [T!fl IT^l H31 150] . 
Recently it has been shown that such a generalization of Gabor frames to general time- 
frequency lattices does not enable us to overcome the Balian-Low theorem [21 E| ■ 
This leads naturally to the question if it is possible to extend the construction of 
Wilson bases to general time-frequency lattices. 

Another motivation for the research presented in this paper has its origin in 
wireless communication. Orthogonal frequency division multiplexing (OFDM) is a 
wireless transmission technology employing a set of transmission functions which is 
usually associated with a rectangular time- frequency lattice |16| . The connection to 
Gabor theory is given by the fact that the collection of transmission pulses in OFDM 
can be interpreted as a Gabor system, see [32] ■ The density of the associated 
rectangular time-frequency lattice can be seen as a measure of the spectral efficiency 
in terms of number of bits transmitted per Hertz per second. The necessary condition 
of linear independence of the transmission functions implies that we are dealing with 
either an undersampled or a critically sampled Gabor system. 

For wireless channels that are time-dispersive (due to multipath) and frequency- 
dispersive (due to the Doppler effect) good time-frequency localization of the trans- 
mission pulses is essential to mitigate the interferences caused by the dispersion of 
the channel |291 135] . The ideal set of transmission pulses should therefore possess 
(i) good time- frequency localization and (ii) maximize the spectral efficiency, i.e., 
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the transmission functions should correspond to a (critically sampled) Gabor basis 
for i 2 (M). As we know, the Balian-Low theorem prohibits these conditions to be 
fulfilled simultaneously. 

Recently it has been shown that in case of time-frequency dispersive channels the 
performance of OFDM can be improved when using general time- frequency lattices, 
in particular hexagonal- type lattices |22j- I n a nutshell, lattices that are adapted to 
the shape of the Wigner distribution of the transmission pulses allow for a better 
"packing" of the time-frequency plane, which in turn can be used to either achieve 
higher data rates or to improve interference robustness of the associated so-called 
Lattice-OFDM system. 

One variation of OFDM (for rectangular lattices) is called Offset-QAM (OQAM) 
OFDM, it corresponds to using a Wilson basis as set of transmission functions [5J. 
OQAM-OFDM achieves maximal spectral efficiency and allows for transmission func- 
tions with good time-frequency localization. As in the case of standard OFDM it 
would be potentially useful for time-frequency dispersive channels to extend OQAM- 
OFDM to general time- frequency lattices in order to improve the robustness of OQAM- 
OFDM against interference even further. Thus we again arrive at the problem of 
constructing Wilson bases for general non-separable time-frequency lattices. 

Yet another motivation comes from filter bank theory, more precisely cosine- 
modulated filter banks [Jj. We know that discrete-time Wilson bases correspond to a 
special class of cosine- modulated filter banks (see [7]). In light of the improvements 
gained by using general time- frequency lattices in OFDM |33| . it would be interesting 
to analyze if the construction of cosine-modulated filter banks can be extended to 
general time-frequency lattices. A positive answer to this question might lead to a 
more efficient encoding of signals and images. 

Since our goal to construct Wilson bases for general time-frequency lattices is 
in part motivated by applied problems and since any numerical implementation of 
Wilson bases is based on a discrete model, our analysis will not only concern L 2 (M.) 
but also comprise the spaces £ 2 (Z) and C L . Furthermore, i 2 {1) is the appropriate 
setting when Wilson bases are utilized as filter banks, since in this case one deals with 
sampled, thus discrete-time, signals. 

1.1. Notation. We assume that the reader is familiar with the theory of Gabor 
frames and refer to [23J for background and details. 

A lattice A in M. d is a discrete subgroup with compact quotient, i.e., there exists 
a matrix A S GL(d, K) such that A = AZ d . The matrix A is called the (non- unique) 
generator matrix for A. The volume of A is vol(A) = | det(A)|. Two lattices, which 
play a crucial role in OFDM design (see jSHDi are the rectangular lattice An and the 
hexagonal lattice Ah- A generator matrix for Ar is given by 



and a generator matrix for Ah is given by 



V2rp y/g rp 



W 2^3 




A H 



where T,F > 0. An easy calculation shows that both lattices A# and Ah have the 
same volume TF. A normal form for matrices, which we will use in the following, 
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is the so-called Hermite normal form 



We say that a matrix A = 



a b 
c d 



is m 



Hermite normal form, if c = 0, a, d > 0, and < b < a. For example both matrices 
Ar and Ah are in Hermite normal form. 

For (x, y) £ R 2 and g £ L 2 (R), let g XjV be denned by 

g x , y (t)=e 2 * iyt 9(t-x). 

We denote by Q(g, A) the system of functions given by 

<?A, M (*)=9(*-A)e 2 ^, (A, p) £ A. 



As usual, the redundancy of G(g, A) is given by vol ^ A ^ ■ 
As in }20| . we define the Schrddinger representation p 



U(L : 



by 



y, z)ff(t) = ^e-^e^ffCt - x). 



Note that 



(1.1) 



Furthermore, we use the following notation from [23j (with slight changes): 



J = 







B h = 



b 
i 



C r .= 



1 
c 1 



J 7 and A denote the Fourier transform, T 1 and v denote the inverse Fourier trans- 
form. The dilation is given by T>tf(t) = \b\? f(bi) and the "chirp" operator is defined 



via Af c f(t) 



f(t), where b, c £ K and / £ L 2 (R). 



Metaplectic transforms will turn out to be a very useful tool in our analysis. For 
the study of the discrete and finite case, we need a result on metaplectic transforms 
from . Since this thesis is not easy to access, we present the result here in the 
slightly weaker version we will use, together with the necessary definitions and nota- 
tions. The result will be stated in the situation of a general locally compact abelian 
group G with dual group 67, group multiplication denoted by +, and action of G on 
G denoted by (x, x) for x £ G and x € 67. The cases we are interested in later on 
are 67 = Z and G = It,. We will usually write a metaplectic transform a in the 

a (3 



matrix notation a — 



7 



£ Hom(6' x 67), which means that a £ Hom(67), 



P £ Hom(G, G),7£ Hom(67, G), and 6 £ Hom(G). Then the adjoint a* is defined by 
{(x,X),<T*(y,n)) = {<r(x,x),(y,n)) for all {x,x),(y,n) £ G x G. Let 77 be defined by 
£ Hom(6' x G,G x G), where the above definition concerning the 



n 



-I? 





la 0" 

matrix notation has to be adapted in an obvious way, and where Ig and Iq denote 
the identity on 67 and G, respectively. Then a is called symplectic, if a*r]a = rj. If 
C £ Hom(6 l x G,G x G), then ip is a second degree character of G x G associated to 

C, if i>((x, x) + (y, t)) = i/>(x, X)i>(y, t) ((%, X), C(y, ^)> for all {x, x), (2/, n) £ G x G. 
Moreover, for (x,x) £ G x G and g £ L 2 (G), let g X:X be defined by 



Jx,X 



(t) = x(t)g(t-x), 



which generalizes the previous definition. 
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We can now state the version of |261 Theorem 1.1.28], which we will employ in 
Sections 01 and 0] 

ad 

6 Hom(G x G) be symplectic and let ip^ be a 



Theorem 1.1. Let a :— 



7 S 



second degree character of G x G associated to 

c ■=<?* 

If U is defined by 






" 


a — 





" 













€ Hom(G x G,G x G). 



Uf(t) 



f(at + Puj)ip- 1 (t,Lu)du;, 



then we have 



(Uf) x , x (t)=^ 1 (x, X )Uf a(XtX) (t), (x, X ) e Gx G. 



2. Wilson bases for general lattices — the continuous case. We first show 
that all lattices in R 2 which are important for applications, such as the rectangular 
lattice, the hexagonal lattice, and lattices whose generator matrix has rational entries, 
possess a uniquely determined matrix in Hermite normal form, which we will call the 
canonical generator matrix. In particular, we characterize exactly those lattices, which 
possess a generator matrix in Hermite normal form. 

We then show that it is possible to construct orthonormal Wilson bases for time- 
frequency lattices A with vol(A) = i, which possess a generator matrix in Hermite 
normal form. In principle Wilson systems can be defined for lattices A with vol(A) ^ 
however so far all known constructions of orthogonal Wilson bases for L 2 (R) are 
strictly tied to lattices with volume \ . In light of this fact throughout this paper a 
Wilson system will always be associated with a time-frequency lattice of volume I. 



Lemma 2.1. Let A be a lattice 



Then the following conditions are equiva- 



lent. 



(i) -P2(A) is discrete, where Pi : M 2 — * K, (x,y) i— ► y. 

(ii) There exists a generator matrix A for A which is in Hermite normal form. 
If one of these conditions is satisfied, the matrix A is uniquely determined. 

Proof. Let 



A' 



be an arbitrary generator matrix for A. 

First we prove that (ii) implies (i). By (ii), there exists a matrix 



.4 



a b 
d 



which is in Hermite normal form and which satisfies AZ? = A. Thus P2(A) 
which yields (i). 

Next we show (i) =>■ (ii). For this, we construct a matrix 



dZ, 
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which satisfies the claimed properties. Without loss of generality we assume that 
d! (if d! = and c' ^ we could change the columns of A'). We begin with the 
following observation. Assume that ^7 is not rational. Since i"2(A) is a non-trivial 
discrete, additive subgroup of K, hence a lattice, there exists s G K\{0} such that 
-p2(A) = sZ. Thus c' = sto and a" = sn for some m,n G Z, a contradiction. This 
implies that the quotient §7 is rational. Setting %, = j, fc, Z G Z with gcd(fc,Z) = 1 
and factoring out 4- , without loss of generality we can assume that A' is of the form 



A' 



with a',b',r G R and c',rf' € Z. Now we proceed as follows. First we set p 
gcd(c', d') > 0. We then obtain 



A' 



± 



rdct(A') 
V 





G A, 



since — , — G Z. Hence we can define a by a 



rdct(A') 
P 



> 0. 



In a second step we compute b and d. Since ^ and — are relative prime, there 

d ' *" — 1 (see [^| Theorem 4.4]). Hence c'm + d'n= p 



exist m, 71 G Z such that ~m + 
and we obtain 



A' 


TO 




r(a'm + b'n) 




n 




rp 



e A. 



Now let fc € Z be chosen in such a way that < r(a'm + b'n) + ka < a and define 
b := ka + r(a'm + 6'n) and d := rp. Without loss of generality we can assume that 
d > 0, since otherwise we just take —to and — n instead of to and n. Then 



r(a!m + b'n) 
rp 



k 



G A 



and |ad| 



rdct(A ) r p _ c Jet(^4_ / ) I . This proves that A generates A and is in 
Hermite normal form. 

At last we prove that the matrix in condition (ii) is uniquely determined. For this, 
assume there exist a, b, d, a', b', d' G K with a, a', d, d! > 0, < b < a and < b' < a' 
such that 



A = AI? = A'Z 2 , 



(2.1) 



where 



A 



a b 
d 



and A' 



a' b r 
dl 



By (|2.1|l . there exist to, n G Z with a' = ma + nb and = nd, which implies that 
a' = ma. Again by 1)2.1(1 . we can find k, I G Z such that 6' = ka + lb and a" = Zd. Since 
vol(A) = |arf| = \a'd'\, we obtain |mZ| = 1. Now a,a',d,d' > implies that a = a', 
d = dl , and 1 = 1. Finally, applying this to b' = ka + lb and using that < b < a and 
< b' < a 1 yields b = b'. Thus we have shown A — A' , which completes the proof. □ 
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In the following we restrict our attention to lattices, which possess a generator 
matrix in Hermite normal form. All results in the situation L 2 (M.) could be derived (in 
the same manner, but with much more technical effort) for general lattices, however 
with little or no practical benefit. 

Definition 2.2. Let A be a lattice inM 2 , which satisfies the conditions of Lemma 
\2.1\ Then the uniquely determined generator matrix A of Lemma \2.1\ is called the 
canonical generator matrix for A. 

Now let A be a lattice with vol(A) ~ i, which possesses a generator matrix in 
Hermite normal form. Using the definition of a canonical generator matrix, we define 
a Wilson system associated with A as follows. 

Definition 2.3. IfQ(g,A) C L 2 (R) is a Gabor system of redundancy 2, and 



A = 



a b 
d 



is the canonical generator matrix for the lattice A, then the associated Wilson system 
W(g, A, L 2 (M)) = {-0m n }mez,n>o consists of the functions 



1p m ,Q = 92ma,0, if n = 0, 



n even. 



V>m,„ = ^/fe " " {g ma +nb,nd + 9rna-nb-nd), if Tl ^ 0, 171 + 
^m,n = -±=e~ nbdn2 {g ma +nb,nd ~ 9ma-nb,-nd), if n 7^ 0, m + n odd. 

If the system W(g, A, L 2 (R)) is an orthonormal basis for L 2 (R) we call it a Wilson 
(orthonormal) basis. 

We will sec that this definition reduces to the usual definition of Wilson systems. 
For this, we consider the rectangular lattice 

r = {(f,n)} m ,„ ez . 

It is an easy calculation to show that the canonical generator matrix for T is 



.4 



I 
1 



Thus, for each g e L 2 (R), the Wilson system W(g, T, L 2 (R)) consists indeed of the 
functions 

^m : o = flVn.Oi if n = 0, 

^m,n = 75 (ffm/2,7i + 9m/2,-n), ^m + n is even, 

^m,n = 7f (ffm/2,n - 9m/2,-n), if TO + n is odd, 

which coincides with the usual definition of Wilson systems, cf. for instance |201 Def- 
inition 8.5.1]. Notice that we will fix the notation T for the remainder. 

We will make use of the following well-known theorem about a Wilson system for 
rectangular lattices to constitute an orthonormal basis (e.g. cf. 3, Theorem 4.1]). 

Theorem 2.4. Suppose that g e L 2 (R) is such that 

(a) g is real-valued and 

(b) {jm/2,n}m.n6Z is a tight Gabor frame for L 2 (M) with frame bound 2. 
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Then the system W(g,T, L 2 (R)) is a Wilson orthonormal basis for L 2 {K) . 

We are now ready to extend the construction of Wilson bases for time-frequency 
lattices which possess a generator matrix in Hermite normal form. 

Theorem 2.5. Let A be a lattice in M. 2 with vol(A) = | and canonical generator 
matrix 

A a b 
L° d 

Define U by 

U :=V 1/d oToAf_ b/d oT-\ (2.2) 

Let g € L 2 (R) be such that 
(i) Ug is real-valued, 

(^) {gma-\-nb,nd }m,nez ^ a tight frame for L 2 (R) with frame bound 2. 
Then the system W{g, A, L 2 (R)) is a Wilson orthonormal basis for L • (M) . 

Proof. We will reduce our claim to Theorem l2.4l bv using a metaplectic transform. 
We define A by 



A 



d -b 
2a 



Then, for all m, n £ Z, we have 

A(ma + nb, nd) = (^rn, nj . (2-3) 

Since A G Sp(2,IR) we can apply \H\ Theorem 4.51] and write 

A = B d (-J)C b / d J . 

By |20l Example 9.4.1] we obtain 

p(x, y, l)g = U- x p{A{x, y), l)(Ug), (2.4) 

with U defined in l|2.2jl . Using (|1 . If) . and (|2.4jl . we compute 

gma+nb,nd = e^ ma+nb > d p{ma + nb, nd, l)g 

= e^^+^'^U- 1 p{A{ma + nb, nd), l)(Ug) 
= eMma+nb^djj-l p ^ m ^ l)(Ug). 

Using (|1 . 1|) again, we obtain 

_ „Tri(ma+nb)nd — iri^mnrj — 1{tt,,\ nibdn 2 tj — 1 /rr \ ( r, r\ 

g ma +nb,nd = e e 2 U {Ug) m / 2 ,n — e U {Ug) m / 2 , n - (2.5J 

Since multiplication by a phase factor and applying a unitary operator to a tight 
frame preserves tightness (and frame bounds), it follows that condition (ii) is equiv- 
alent to {{Ug) m /2,n}m,n<EZ being a tight frame with frame bound 2. We need not 
deal with condition (i), since this states already that Ug is real- valued. Applying 
Theorem 12 .41 yields that the Wilson system W(Ug, T, L 2 (K.)) is an orthonormal basis. 
Using now the metaplectic transform, i.e., I|2.5|l . and the fact that U is a unitary 
operator, finishes the proof. □ 
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We will conclude this section by providing an example, which shows how we can 
compute a Wilson basis with excellent time-frequency localization for a special lattice, 
but this calculation can also be done for an arbitrary lattice. 

Here we consider the hexagonal lattice Ah with generator matrix 



A H = 



%/2 %M ' 

» « 



This lattice was also used in [331 - Observe that since we have < < -^|, the 
matrix Ah is already the canonical generator matrix for Ah- First we define the 
function h € L 2 (R) by 



h(x) = {2v)±e 



By (201 Theorem 7.5.3], the set {h m / 2 ,n}m,nez is a frame for L 2 (R). Let S denote its 
frame operator and consider the function ip G L 2 (R) given by 



'•P 



= V2To S~*h. (2.6) 



Using [3] Theorem 4.6], this function coincides with the function considered in |1 II 
Section 4]. There it was shown that <p satisfies conditions (a) and (b) of Theorem l2.4l 
and hence yields a Wilson basis in the sense of Theorcm l2.4l Moreover, ip has expo- 
nential decay in time and frequency. To obtain a generating function for a Wilson 
basis with respect to Ah in the sense of Theorem 12. 51 first observe that, by the proof 
of Theorem 12 .51 we only need to compute the function g = U^ 1 ^, where U is defined 
in (|2.2|l . Then g automatically satisfies conditions (i) and (ii) of Theorem 12.51 and 
hence the system W(g, Ah, L 2 (M.)) is a Wilson orthonormal basis by Theorem 12.51 
Thus we define g € L 2 (R) by 

g = FoN i o T~ x oVj^tp = \/2ToNj_ oV^o S~^h. 

Let us mention that the function g has exponential decay in time and frequency. 
Thus we obtain a Wilson basis with respect to the lattice A with very good time- 
frequency localization. As already mentioned above this procedure can be applied 
to an arbitrary lattice, hence we obtain a Wilson basis with excellent time-frequency 
localization for any lattice. 

3. Wilson bases for general lattices — the discrete case. In this section we 
analyze the construction of Wilson bases for general time-frequency lattices for func- 
tions defined on £ 2 (Z). The reasons for considering the setting £ 2 (Z) are on that the 
one hand several applications such as filter bank design in digital signal processing deal 
directly with a discrete setting [7j, and on the other hand, even those problems that 
arise in the "continuous" setting of L 2 (R) require a discrete model for their numerical 
treatment. Thus, with these practical aspects in mind, throughout this section we 
naturally consider only lattices whose generator matrices have rational entries, since 
any implementation is intrinsically restricted to such "rationally" generated lattices. 
Another natural setting for numerical implementations is of course C L (which can be 
identified with the space of L-periodic sequences). We will analyze that case in the 
next section. 
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Before we proceed we define Gabor systems and Wilson systems on £ 2 (Z) for 
general time-frequency lattices A with vol(A) = |. First we prove that each lattice 
possesses a generator matrix of some particular form. 

Lemma 3.1. Let A be a lattice in Z x K with generator matrix A given by 



A = 



with a, b E Z, c, d E Q, and det(A) = — , 



and denote c — = |j wii/i gcd(p, g) = gcd(p',q') — 1, p,q,p',q' E Z. TTien A 
possesses a uniquely determined generator matrix of the form 



A' 



2. b' 

2 " 



(3.1) 



where N = 



and b' E Z, < 6' < f 



gcd(pg',p'g) " ^ " — " ^2 

Proof. We assume that c ^ 0, otherwise (|3.1(l is automatically satisfied. We first 
show that A can be written as 



.4 



(3.2) 



a b 

r s 
. N N _ 

with integers r, s, TV, such that gcd(r, s) = 1. 

Let c = £,d = 2- with p,p',q,q' E Z, denote TV' := gg',c := := p'g and 

write z = gcd(c, d). Since a,b,c,dE Z and since 

1 N' 
vol(A) = — => ad — be = — , 



it follows that y£Z. A necessary and sufficient condition for the equation ad—bc = 
2- to have an integer solution in a and b is that gcd(c, d)\2-, see j^J Theorem 8.1], 
hence z\2-. Denote z' :— S-,r := -,s := -. Then c = ^r.d = k^t with z' El and 

1 2 2.2 ' z ■ z 2z' ' 2z 

gcd(r, s) = 1. By a proper choice of the signs of c and d we can always assume that z' 
is positive. By writing N := 2z' E Z we see that ^4 can indeed be written as in i|3.2|) . 
Now, assuming that A is of the form (|3.2(l . we compute 



a 6 




' TVd" 




~N(ad - be) 




"AT 
2 


c d 




-Nc 













Since gcd(r, s) = 1 there exist integers m, n with mr ■ 
we denote b' = am + bn and obtain 



1. For such a pair (m, n) 



a 6 




TO 




a?n + bn 




~b'~ 


C G? 




n 








1 

.N. 



If 6' < or b' > y, we substitute the vector obtained by 



~b'~ 

1 


+ k 


"JV" 

2 





'6' + fcf " 
1 z 


. JV - 






N 



where k E Z is chosen in such a way that < b' + k^ < 2. Consequently the matrix 



(3.3) 
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generates the lattice A. Finally, since y and b' are integers and < b' < y, the 
generator matrix in (|3.3[1 is indeed of the form l|3.1|l . 

The fact that this is a unique representation follows immediately from the condi- 
tion 0< b' < f . □ 

Definition 3.2. Let A be a lattice in Z x R. TTien i/ie uniquely determined 
matrix A' of Lemma \S.l\ is called the canonical generator matrix for A. 

In the following we will regard such a lattice as a lattice in Z x T by considering 
A = {f m + bn,e 2nl %} mez,n=o,...,JV-i- This is a very natural approach, since, for all 
k G Z, we have 

(f m + bn, jjn + k) = (f (to - 2bk) + &(n + kN), j?(n + kN)), 

hence the lattice A' J? is invariant under adding integers to the second component. 
Moreover, it is sufficient to restrict to the index set Z x {0, . . . , N — 1}, since, for all 
< ri < N and k € Z, 

(fm + b(n' + kN),(n' + kN) mod iV) = (f (m + 2Mc) + 6n',n'), 

which implies 

{(^m + bn,n mod iV)}m,nez = {(fm + 6n,n)} me z,„ = o,...,jv-i 
in the sense of sets. 

Using the definition of canonical generator matrices we can now define Gabor 
systems for £ 2 (Z). 

Definition 3.3. Let A be a lattice in Z X T with canonical generator matrix A 
given by 



and let g S ^(Z). Then the associated Gabor system {g m N_ +nb n i_} m ei,,n=o,....N-i is 
given by 

9m K +nb<njr (l) = g(l - (mf + nb))e Mln / N , I G Z. 



Now we first give the definition of a Wilson basis associated with a lattice with 
diagonal canonical generator matrix, i.e., with 6 = (in this special case the definition 
coincides with the one given in 0). 

Definition 3.4. Let A be a lattice in Z X T with canonical generator matrix A 
given by 

f o 

and let g £ ^(Z). Then the Wilson system W(g, A,^a(Z)) = {V'm,n} TO gz n=o ^ * s 
given by 



and for togZ, n = 1, . . . , y — 1, 



A 

m,n 
A 

m.n 



V ,.JL +3 m f , 



N „ 1 



z/ m G Z, n = 0, y , 

if m-\- n even, 
if m + n odd. 
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The Zak transform, which can be defined for any locally compact abelian group 
(cf. HHDi will be employed to prove equivalent conditions for the Wilson system to 
form an orthonormal basis. In particular, we need the Zak transform on T with 
respect to the uniform lattice K = {e 27 ™«" : k = 0, . . . , ^ — 1} in T, which is defined 
on the set of square-integrable functions on {e 27 ™* : t S [0, j^)} x {0, . . . , y — 1} by 



Zf(e 2mt ,y) = J2 f(e 2 ^ t+ ^)e 2m ^y. 

k=0 

The proof of the following proposition is inspired by the proof of [111 Proposition 

5.2]. 

Proposition 3.5. Let g € £2(^1) be such that g is real-valued and consider the 
lattice A with canonical generator matrix given by 

jr. 

Then the following conditions are equivalent. 

(1) idm^- n i}mez,n=o....,JV-i is a tight frame for £2^) with frame bound 2. 

(ii) We have \Zg(e 27rlt , y)\ 2 + \Zg{e 27rz( - t+ ir\y)\ 2 = N a.e.. 

(iii) For all j€{0,..., N-l}, we have X^J 1 g(e 2 ^ t+ ^)g{e 2 ^ t+ ^) = NS j>0 
a.e.. 

(iv) W(g, A, £2^)) is an orthonormal basis for li^K). 

Proof. Throughout this proof we choose the normalized Haar measure on T, i.e., 
m(E) = J Q l£;(e 27 ™*)ctt for all measurable E C T and the counting measure on its 

dual group T = 1. This choice ensures that the Plancherel formula for T holds. 

Since we will mainly work in the Fourier domain, we first need to compute the 
Fourier transform of the elements of the Gabor system for our following calculations: 

{e 2mt ) =J2e 2 ™% l g(e 2 ™«- m ^)e- 2 ™ H 



I 1 \mn * ( 2-nit \ 



(i) (ii): Since the Fourier transform is a unitary operator, the Gabor system 
{<7 m w j_„}mez,n=o....,JV-i is a tight frame for £2^) with frame bound 2 if and only 
if the Gabor system {g n ^ m » }jn,ez,n=o,...,JV-i is a tight frame for L 2 (T) with frame 
bound 2. Then we write this set as the disjoint union {g n 2_ m "} m gzn=o —-1 U 
{(T_±g) n 2_ m N_ } mgZ „ =0 jv„ 1 =: Gi UG2, and let Si denote the frame operator for 
G h i = l,l' 2 ' ' 

First we compute the frame operator Si. For all / € L 2 (T), using the Poisson 
summation formula O Theorem 4.42] applied to H = {e 2lTlk ~ : k — 0, . . . , ^ — 1}; 
we obtain 



N 



1 



Sl f(e 2 ^) = £ E >*nW 

mGZ n=0 
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N , 

pi 

= EE / f(e 27Tls )g{e 2m ( s - n 7T)) e - 2mm % s ds 3(e 2 " (t ""^ ) )e 27r4m ^ t 



-1 



E 



2 (/T 4 j)(mf)e 2 "^' 



n— Lm 



JV i JV 



-1 4f-l 



n=0 fc=0 

Applying now the Zak transform yields 



Z(S 1 /)(e 2 - t ,y)= E ^ E /(e 2 " (t+ ^ ) )g( e 2 "( t + 2ii± ^))5(e 2 " (t+2i ^ ) ) e 2T ^ ! ' 

;=o n,fc=0 



= 2, V /(e 2 "( t + 2i ^))g(e 2 "( t +#))g(e 2 "( t+2i ^^)e 2 ^ 2 " " + " ' « 

(,n,fc=0 



= ^E .9(^ (t+#) ) 

_1 



g{e 2m{t+ ^^)e~ 2 ^y 

k=0 



n=0 



= j } Z{g) (e^ , y)Z(g) (e 2 ^ , y)Z(f) (e 2 ^ , y). 

To compute the Zak transform of 52, we can use the previous calculation, which 
yields 

Z(S 2 f)(e M \y) = ^Z(.f){e 2 ^ iy )\Z(T^g){e 2 ^,y)\ 2 
= j r Z(f)(e 2 ™\y)\Z(g)(e 2 ^ t+ ^,y)\ 2 , 

since 



Z(T_^g)(e™,y) = ]T 5 ( e 2 -('+M)) e 2 ™ty = Z{g){e 2 ^^\y). 

fe=0 

Hence (i) is equivalent to 

2Z(f)(e 2 " it , y) = Z((5i + S 2 )f)(e 2 ™\ y) 

= ±Z{f){e 2 ^\y) [\Z(g)(e 2 ^,y)\ 2 + \Z(g)(e 2 ^ t+ ^,y)\ 2 ] a.e., 

which holds if and only if (ii) is satisfied. 

(ii) (iii) : The following properties of the Zak transform will be exploit several 
times. The reconstruction formula 



n 1 



n 1 



Zf{e 2 ™\y) = f(e 2m(t+¥) ) E e 2 ^ = f f{e 2mt ) 



fc=0 



j/=0 
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holds a.e., since eJLq 1 e 2 ™^ ^ if and only if k = by [HI Lemma 23.29] and, if 
k = 0, then X^o* e 2jrl ^~ y — Moreover, we will use that 



Zg(e 2 ™( t+ %\y) = fie^+^ly 2 ^ = e' 2 ^ Zg(e 2mt , y). 

k=0 

The idea is to write the equation in (iii) in terms of the Zak transform. Using the 
fact that g is real-valued, we compute 



JV-l 



2™(*+^))£( e 2™(t+^£) 



1=0 



JV-lf-1 



N 2 



E E Zg{e 2 ™^\x) J2 Zg(e^+^),y) 



1=0 x=0 

o - 1 - 



iV 2 



E [^(e 2 " (t+ -),x)Z5(e 2 -(*+^),y) 



k,x,y— 



+Zg(e 2m ( t+ ^),x)Z5(e 2 -( t + i ^ ±i ),y) 



iV 2 ^ 

x,y— 



fc=0 



.2^22. j, 



Z5(e 2 " t , a; )Zg(e 2 - t ,y) 



+Z3(e 2 " (t+ ^ ) ,.x)Z5(e 2 "( t +^),y) 
f- 1 

= Jf E [l^ 2 "*^)! 2 + \Z 9 {e 2 ^^\x)\ 2 

x=0 

= | [|^(e 2 ^, -)| 2 + |Zg(e 2 ^ t+ -), -)| 2 ] " (j), 

where the inverse Fourier transform is taken in Z at . This shows that (iii) is equivalent 

T 

to 



N 



\Zg(e 2 ^, -)| 2 + \Zg(e 2 ^ t+ ^\ -)\ 2 } (j) = N6 jfi . (3.4) 



If (ii) holds, then 
2 



N 



\Zg(e 2 ™\ .)| 2 + \Zg(e 2 ^+i\ -)\ 2 (j) = ±N E = A^.o, 



which is 1|3.4[) . On the other hand, the inverse Fourier transform is injective. This 
proves that i|3.4[l holds if and only if (ii) is true and thus (ii) (iii). 

(iii) (iv): First we remark that W(g, A, liffij) is an orthonormal basis if and 
only if the set 



* :={T nN f m :m = l,...,N, n e Z}, 
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where 

fi(x) = g(x), 
Jn{x) =9 .m{x), 

= ^L-(g k§>jr + (-l) k+l g kf ^x), I = 1, . . . , f - 1, k = 0, 1 

is an orthonormal basis, since these elements differ from the elements in W(g, A, £ 2 (Z)) 
only by factors of absolute value 1. Next notice that to prove (iv) it is sufficient and 
necessary that 

||T„iv/ m ||2 = 1, m=l,...,I,neZ (3.5) 

and 

L 

/) /, {hi, T nN f m ) (T nN f m , h 2 ) = (hi, h 2 ) , for all hi, h 2 S £ 2 {Z). (3.6) 

m=l ngZ 

We start by dealing with (|3.5|l . Using the Plancherel theorem, we compute 
1 = || Willi - Ml = I 1 g(e 2 « a )g{e 2 ^)dt, 



1 = ||W;v||I = ||.g« || 2 - f 1 g(e Mt )g(e 2 ^)dt 



and, for m = 2, . . . , N — 1, 

1 = ||7nJv/m||2 

= ll^(%- fcf + (-l) fe+i .9- 1 L^ f )||^ 

1 /- 1 







1 . 







|g( e 2 "(*-ir))|2 + |^ (e 2 7 ri(t+^) ) |2 + (_ 1 )fe+!^( e 2 W i(t-l)^ (e 2 7 ri( t +J r ) ) 



+ (-l) fe+ 'g(e 27rl ( t -^))g(e 2 "( t +^)) 
Since g is real-valued, we can continue the last computation and obtain that 

1 = \\T nN f m \\ 2 2 = \\g\\l + (~l) k+l f 1 g(e 2 ^)g(e 2 ^ t+ ^)dt. 

Jo 

Combining the above computations we have proven that 1|3.5[) holds if and only if 

f 1 g(e 2mt )g(e 2 ^ t+ ^)dt = 5 jfi for all j € {0, . . . , N - 1}. (3.7) 
Jo 

Now we turn to the study of condition Ij3.6|l . Using the Plancherel formula and 
the Poisson summation formula Theorem 4.42] applied to H = {e 27 ™"^ : k = 
0, . . . , N — 1}, we obtain 

N 

^ ^ (hi,T nN f m ) (T nN f m , h 2 ) 

m=l n£Z 
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N 



m—1 n£ 
N „i 



E E (h\,T nN f m j (T nN f m , h 2 

rn—1 n£Z 

E E / (hu^™) / (£M(e 2 ™) 

E(^H Ar ") e2 " tAr " 



2 ™ sNn ds e 2mtNn dt 



.nSZ 
AT-1 



r=0 



= E / (hiU(e^) 

= E / (^/™)(e 2m *)^E(/™^)(e 2 " (t+ ^))^, 
which equals (/ii, /12) if and only if 



dt 



N 



J2f"^ 2mt )fm(e 2m(t+ ^)^N6 rfi foraUre{0,...,JV-l}. (3.6 



Setting L := {— -j + !,■•■, — 1, 1, . . . , ^ — 1}, we compute 



m—1 



j(^(^*)) + j f0 (e«)j f0 (e^*))t^ E^,-fef( e27rit ) 



=1 fc=0 



1 f_1 1 

g(e 27 "*)g(e 2 ™( t+ ^) + g(e 2 ^-f )).9(e 27r4 ( t+ ^-f >) + - E E 



1=1 k=0 



+ (-l) k+l g(e 2 ^ t+ ^)}[g(e 2 ^ t -^ + ^) + (-l) k+l g( e 2m{t+ ^+^)] 
= g(e 2m *)g(e 2 "( t+ ^)) + ^e^-T^e 2 ^*-^) 

+ 1 ^ ^(_l)fer [g( e 2^(t-17))^( e 2«(t-i7 + #)) + g( e 2^(t+17))g( e 2«(t+17 + ^)) 

^ ;=i k=o 



- g{e 2mt )g(e 27Tl( - t+ 

+ E(- 1 )'9( e2 " (t+i) )9( e2 " (t_i+ ^ ) ) 



) )+Eff( e 



2"(*+w))g( e 2"(t+i + w)) 



i+(-ir 



ZGL 



l + (-l) r+1 



+ ^ (e 2 m(t -f)^ (e 2«( t+ ^-f) 



If r is even, i.e., r = 2j, we obtain 



N-l 



E /™(e 2 ^)/ m (e 2 " ( * + ^)) = E .g(e 2 ™ (t+ ^)g(e 2 "( t+ w+#)) 7 



m—1 



Z=0 
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and if r is odd, i.e., r = 2j + 1, we obtain 

jv 



m=l 

= g(e 2mt )g(e 2 ^ t+ ^) + ^^'-^(^-(t+^f )) 
+ J2(-l) l g(e 2mt - t+ Try)g(e 27r ^ t ~7r + ^) 



This shows that (|3.6I) holds if and only if 



N-l 



J2 g{e 2 ^ t+ ^)g{e 2 ^ t+ ^) = NS j}0 . 



1=0 



Moreover, this equation implies equation (|3.7|l . since 



l(e 2mt )g(e 2 ^ t+ ^)dt = 



1=0 



2 ^(*+w))g( e 2 "( t +TT i )W^ = 



)dt = / N5 3fi dt = S jfi - 



This shows (iii) O- (iv), and hence the theorem is proven. □ 

Now we will study the case of a general time-frequency lattice. 

Proposition 3.6. Let A be a lattice in Z x T with vol(A) = \ with canonical 

generator matrix A given by 



A 



£ b 
2 " 

A7 



Let g £ ^(Z), let mo,rto G Z &e chosen such that Y m o + ^"o = gcd(-jifr) — : c ; arl ^ 
Ze£ [/ &e defined on ^(Z) by 



Then 



Uf(k) = f(k)e^ k2 . 



,,N +nbrn ^(l) = C(m,n)U(U 1 g) m N +nb ^ m ^ n ^ +n Al), leZ, 



where C(m,n) = e ™^r(™f+nb)\ 

Proof. Let a £ Hom(Z x T) be defined by 



h 

_ JUL U 
cN 1T 



It is easy to check that a is symplectic on Z x T. In order to apply Theorem ll.il we 
need to compute a second degree character of Z x T associated to C, = a* k^u — Kq, 
where kq is defined by 



K 



1 1 

h 



G Hom(Z x T, T x Z) 



First we note that a*(z,n) = (ze 2m ^ since 

((m,t),a*(z,n)) = (a(m,t),{z,n)) = z m (e 2 ^-^ m t) n = <J(m, t), {ze 2m{ -^ )n , n)\ 
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Thus 

C(n, *) = {a*K a - K )(n, z) = a* K (n, e 2 ^^ )n z) - (l,n) = (e 2 ™ ( ~^ )n , 0). 
The map 

fZxT^T, ip(m, t) = e - 7Tt ^ m2 
is a second degree character associated to £ as the following calculation shows: 

ip(m,t)ip(n, z) ((m,i),((n,z)) = e -™w™V™^ n2 e 27ri( -^ mn 

= ^((m,i) + (ra,z)). 

Next notice that 

CT (mf + n&, 4) = (mf + n&, -to§* - + n±). (3.9) 

Now we can apply Theorcm ll.il which proves the claim. □ 

Next we define a Wilson basis associated with a lattice with arbitrary canonical 
generator matrix. For this, the following mapping will turn out to be very useful. 

Lemma 3.7. Let y,6 € Z with < b < y, and let mo,rto € Z be chosen such 
that jmo + brio = g°d( : y,6) =: c. Further let d := lcm(^r,6). TTien t/ie mapping 
if : Z 2 — > Z 2 defined by 

(m n) = I (m, n) : 6 = 0, 

^ ' ' [ (mm — ^n, ran + |n) : 7^ 

is bijective and, for all to € Z, we have 

{(m, n mod 2c) : (m, n) G p -1 (Z x {0, . . . , N - 1})} = {m} x {0, . . . , 2c - 1} 

with 

\{n : (to, n) E ^ _1 (Z x {0, . . . , N - 1})}| = 2c. 



Proof. We only need to study the case b ^ 0. For this, let (m,n), (m',n') G Z 2 
be such that <p(m, n) = cp(m' ', n'). Then 

§{mm - jj-n) + b(mn + ±n) = f (m'm a ~ ^n) + b(m'n + %n'), 

which holds if and only if 

and hence to = to'. This implies 

(-^n,fn) = (-fn' ) fn') ) 

which yields n = n' . This proves that tp is injective. 

To show that ip is surjective, let (k,l) 6 Z 2 and consider M := yfe + W. It is 
well-known that there exists some to 6 Z with M = toc. Furthermore, we have 

{(p, g) G Z 2 : -yP + = mc} = {(totoq — fy-n, mno + ^n) : n € Z}, 
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since + bq= %P' + W if an d only if y (p — p') = b(q' — q). This yields the existence 
of some n£Z with 

y(m, n) = (mm — 2^n, mn + f ft) == (k, I)- 

Secondly, we will prove the second part of the lemma. First observe that ra,neZ 
satisfy 

<p(m,n) £ Zx {0,...,N- 1} (3.10) 

if and only if they satisfy 

-gnom < n < -±n Q m + |(/V - 1) = 2c - |n m - |. 

Hence, for each fixed to € Z, the set of n S Z such that l|3.10[) is satisfied equals 

Sm ■= {l-^nom], . . . , [2c- ±n m - jjj}. 

To finish the proof we claim that 

|5' m |=2c for all to eZ. (3.11) 

For this, fix to £ Z and let k £ Z and Z £ {0, — 1} be such that —n^m = fc| + /. 

Then we obtain [~— ^n TO,] = if i = and otherwise [— ^n rn] = fc + 1. Moreover, 
we have 



L 2 c -|noTO-|j 



2c+fc+ ^± 



which equals 2c + — 1 if / = and otherwise 2c + fc. Thus the second part of the 
lemma is proven. □ 

Note that the following definition reduces to Definition ^. 4l in the case of a diagonal 
canonical generator matrix. 

Definition 3.8. Let A be a lattice in Z x T with canonical generator matrix A 
given by 



A 



& b 

2 " 



Let g £ ^(Z), and let ip be defined as in Lemma \3. 7[ Then the Wilson system 
W(g, A,£ 2 (Z)) = {Vw}mez,n=o,...,# is 9™en by 



^m,n 9ip 1 (2m,n)-^--|-cp2(2m,n)fc,^2(2m,n)-^ 



if to £ Z, n = 0, -y , 



and for raeZ, n = 1, . 



iV 



_), if m + n even, 



i(m,n)-^r-+^2(fn,™)fr,</>2 (m.n) - 
VVi,n 75 (^¥'l(m,n)-j-+i/32(m.n)f),¥'2(m,n)-^ 

-ff^^m.-njf +1/32 ( m _„)b, V2 ( m ,_„)^), ifm + n odd. 
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The following theorem gives an equivalent condition for a Wilson system with 
respect to an arbitrary time-frequency lattice to form an orthonormal basis in terms 
of a frame condition for the associated Gabor system. 

Theorem 3.9. Let A be a lattice in Z x T with canonical generator matrix A 
given by 



A 







Let g € ^2(Z) be such that U^ 1 g is real-valued, let M := 2c, and let U and ip be defined 
as in Provosition \S. 61 and Lemma \3. 7\ respectively. Then the following conditions are 
equivalent. 

(i) {Sm^+ra6,n 1 }meZ,n=o N-i * s a tight frame for £2(2) with frame bound 2. 

(ii) {((/ _1 ff) m M n i}ma,n=o,...i/-i is a tight frame for £2 (Z) with frame bound 
2. 

(iii) W(f/^9,fZx Aj{0,...,M -1},£ 2 (Z)) is an orthonormal basis for £ 2 {Z) . 

(iv) W(<7, A, £ 2 (Z)) * s ffl?1 orthonormal basis for £ 2 (Z) . 

Proof. Let er be defined as in the proof of Proposition 13.61 and let ip — ((pi,ip 2 ). 
Then we compute 

a((pi (m,n)% + <£ 2 (m,n)&, <£ 2 (ro,n)^) 

= (<pi(m,n)% + •P2{m,n)b,-ipi{m,n)'^ - Lp 2 (m,n)^$ +ip 2 (m,n)±) 
= m(f m + bno) + n(-%% + f &),m(-f m - ^n + ±n ) 

1 /2d no _ d bn , d 

W 2c b cN ^~ b N >> 



K™o(f - ^(fwo + bno)) + 

, M 



where in the last step we used cd = ^-b. Using Lemma [3.71 the equivalence of (i) 
and (ii) now follows immediately from Proposition 13.61 and Equation (|3.9|l . since U is 
unitary and |C(m, n)| = 1. Proposition ^ . 51 proves (ii) -£4> (iii). Therefore it remains to 
prove the equivalence of (iii) and (iv). For this, we will use the following implication 
of Proposition ^. 61 

U(U 1 g) m M n i_ = U(U 1 g) cr ( ¥ , 1 ( m ,„)« +ip2 ( m: „)b !l/ , 2 ( m! „)^) 

= C(<p(m,n)) 5 v , 1 (m,n)^+*'2(m,n)6,vsa(m,n)^- 

Further notice that C((p(m, n))^ 1 does not depend on the sign of n, since 

C(ip(m, n))^ 1 = e _7ri ^(vi("i:")-f + 'P2("i,n)6) 2 _ e —rri^-m 2 (m f +n„b) 2 _ ^--xl^m 2 c 2 _ 

Using now the definition of a Wilson basis, the fact that U is a unitary operator, and 
the fact that \C((p(m,n))\ = 1 yields the result. □ 

4. Wilson bases for general lattices — the finite case. The space C L has 
several advantages over £ 2 (Z) when constructing numerical methods for practical time- 
frequency analysis, which often allow a further acceleration of numerical algorithms, 
e.g., see [Tim]. 

Before defining Gabor systems and Wilson systems for C L for general time- 
frequency lattices, we first prove that each such lattice does not only possess a uniquely 
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determined generator matrix in Hermite normal form (which was already proved in 
|22|L but moreover in our situation this matrix attains a special form. 

Lemma 4.1. Let A be a lattice in Z^ x Z^ with generator matrix A given by 



a b 
c d 



with a, b, c, d £ N, and det(A) 



L 
2' 



and denote p = gcd(c, d) i/c ^ and p = d if c = 0. Then A possesses a uniquely 
determined generator matrix of the form 



A' 



±- b' 

2p 

p 



(4.1) 



where p = gcd(c, d) and < b' < 

Proof. Let p = gcd(c, d) with c 
ad — be = (ar — bq)p = ^. Si"" " 



•mce - = r and — 

p p 



qp : d = rp and note that p\^ since det(A) 
— q, we have 



r 




■ L - 

2p 


-Q_ 








(4.2) 



Furthermore, we claim that there exists a z £ 1 with < z < ^ such that the point 



belongs to A. This can be seen as follows: The condition 
to the existence ofra,n£Z such that 



€ A is equivalent 



a 


b 




m 




z 


c 


d 




n 




P 



(4.3) 



Consider the equation cm + dn — p and substitute c = qp, d = rp, then qpm + rpn = p, 
hence qm + rn = 1. Since r and q are relative prime, there exist m, n £ Z such that 
qm + rn = 1 (see e.g. [21 Theorem 4.4]). Thus (14 holds for z £ Z, but we still 
have to show that it holds under the condition < z < -k-. We can write z = b' + k4 L 
with < b' < tt- and k £ Z. Hence 

— 2p 



z 




V 


+ k 


■ L - 

2p 


p 




p 








Since 



- L - 






2p 


= A 


r 







_-q_ 



by JO), it follows that 



£ A. Consequently the matrix 



A' = 



2p 

p 



(which satisfies det(A') = -|) generates A. 

The fact that this matrix is uniquely determined is an immediate consequence 
from the condition < b' < □ 

— 2p 

Definition 4.2. Let A be a lattice in 7Ll x Z^. Then the uniquely determined 
matrix A' of Lemma \j^.l\ is called the canonical generator matrix for A. 



WILSON BASES FOR GENERAL TIME-FREQUENCY LATTICES 



21 



Using the notion of a canonical generator matrix, we first give the definition of a 
Gabor system. 

Definition 4.3. Let A be a lattice in Z^ x 7L>l with canonical generator matrix 
A given by 



A 



2p 

p 



Set M = 2p, N = - and let g be some L-periodic function on Z. Then the associated 
Gabor system is given by {g ma +nb,nd}m=o,...,M-hn=0,...,N-i, where 

9ma+nb,nd( l ) = d( l ~ ( ma + nb) )e 2mlnd ' L , I = 0, . . . , L - 1. 



Next we define a Wilson basis associated with a lattice with diagonal canonical 
generator matrix in the following way: 

Definition 4.4. Let A be a lattice in !Ll x Zz, with canonical generator matrix 
A given by 



A 



2p u 

p 



and let g be some L-periodic function on Z. Then the Wilson system W(g, A, C L 
{■*Pm,n}(m,n)ei, where 1 = {0, ...,p- 1} X {0, ^}U{0, ...,2p- 1} X {1, 



is given by 



2p 



yra— ,np ' 



if m = 0, 



,p- 1, n = 0, ^, 



and /or m = 0, . . . , 1p — 1, n = 1, . . . , — 1, 



A 
A 



-^(g m ^, np + g m L_^ np ), ifm + neven, 

-^(g m L-, np - g m L_ : - np ), ifm + n odd. 



Also in the finite case we will employ the Zak transform. This time we will use 
the Zak transform on the group Z^ with respect to the uniform lattice K = {2pk : 
k = 0, . . . , ^ — 1} in Z^, which is defined on the set of square-integrable functions 
on the set {0, . . . , 2p - 1} x {0, . . . , ^ - 1} by 

Zf(x,y)= f(x + 2pk)e 2 ^y, 

k=0 

where we associate Z^ with {0, . . . , L — 1}. 

The following proposition is the analog to Proposition 13 . 51 for the space C L . 

Proposition 4.5. Let g be some L-periodic function on Z such that g is real- 
valued and consider the lattice A with canonical generator matrix given by 

2p 

p 

Then the following conditions are equivalent. 
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(i) {fl , ro ^,np}m=0,...,2p-i,n=0 J ...,f-i is a tight frame for C L with frame bound 2. 

(ii) We have \Zg(x,y)\ 2 + \Zg{x+p,y)\ 2 = ± a.e.. 

— — 1 

(iii) For a// j = 0,..., ^-1 and y S Z L , we have Y,Lo 9{v + l P)g{y + l P + ^iP) = 

(iv) W(g,A,C L ) is an orthonormal basis for C L . 

Proof. The proof, while lengthy, is very similar to the proof of Proposition ^. 51 In 
fact, with obvious adaptations, such as using the normalized Haar measure on Z^, i.e., 
m{E) = j;J2xez L 1e( x ) for all £ C 2^, and replacing Zak transforms and Fourier 
transforms by their corresponding finite counterparts, the proof carries over almost 
line by line. We therefore leave this part to the reader. □ 

Now we will turn our attention to general time-frequency lattices. Here the 
situation is slightly more involved compared to i 2 (Z) . 

Let A be a lattice in Zi x Zl with canonical generator matrix A given by 



.4 



2p 

p 



Then we choose a, [3, to.q, rig € Z such that 



f3pn 



(4.4) 



attains its minimal positive value. Assume that there exists a choice of a, /?, too, no 
such that (^too + bn )(pn ) < and (a^)(ab + (3p) > 0. In the other cases we have 
to change the signs of the later defined 7 and S accordingly. In the following we will 
restrict to the case where |a| = 1. For the remainder of this section let a,/3, mo, no 
be defined in this way. Now we regard a and (3 as elements of Z/,. For the sake of 
brevity we set 



c := gcd(a j 1 , ab + (3p), d := lcm(a ab + (3p) 



and 



s := gcd(^TO + bn a ,pn ), t := lcm(^m + bn ,pn ). 
The minimality condition for (|4.4|) shows that 

a 3p m ° a ^ n o + fiP n o = c = s. 



(4.5) 



We further define 7, S £ Zl by 



7 ^ m o + bn o 



and 5 := — - 



pn 



and a € Hom(Zi x Z^) by 



a 
7 S 



Proposition 4.6. Let A be a lattice in Z^ x wii/i canonical generator matrix 
A given by 
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Let a be defined as in the preceding paragraph, and let U on the space of L -periodic 
functions on 1 be defined by 

Uf(k)= J2 /(afc + /30e^ l(Q7fe2+ ^ 2)(L+1)/L e- 2 ^ 7fe ' /L . 
zez L 

Then 

g m L. +nbtnp (l) = C(m,n)U(U- 1 g) <r{m L. +tlitnp) {l), 

Where C(m,n) = e -^K^ t -m, 2 +PSn 2 )(L+i)/L e -2^i3y m n/L ^ 

Before moving on to the proof of this statement we point out that the operator U 
in Proposition ^. 61 is no longer a simple chirp operator as for the case i 2 (Z), cf. Propo- 
sition l3~BI This difference and the different form of a necessitates a somewhat different 
proof for the case C L . 

Proof. In order to apply Theorcm ll.il we need to check whether a — (01,02) is 
symplectic. For this, we have to show that, for all (x, y), (x' , y') G x Z^, 



e -2iriij 2 (x,y)a 1 (x' ,y') / L ^2ni<j 2 (x' ' ,y')a 1 (x,y) / "L _ ^-2-nix' y / L ^2-nixy' / L 

We have 

e -2Tri<j 2 (x,y)a 1 (x' ,y') / L & 2nia 2 (x' ,y')a t (x,y) / L 
_ e -2TTi(fx+Sy)(ax'+l3y')/L e 2Tri(jx'+Sy')(ax+f}y)/L 
_ e 2TTi(aS-f3i)(xy' — x' y) / L 



(4.6) 



and 



r a Oct §t 

ao — P7 = j 



P n o j^m Q + bn Q 

-t 



Jl , , u :(a(|»no + bn ) + (3pn ) 



By E3, 



aS — /?7 = — j — = 1, 

(j^ m o + bn )(pn ) 



since (^mo + bno)(pno) < 0. This proves that (|4.6|l is satisfied, which shows that o is 
indeed symplectic. Moreover, a short computation analogous to the one in the proof 
of Proposition 13 . 61 shows that 

£(fc, I) = (o* kq<j — Ko)(k, == (ct^fk + (3jl, f3jk + f38l). 

Now it is easy to check (compare also j^HI Example 1.1.34 (hi)]) that 

i> : Z| -> T, rp{k, I) = e ^(^ + [38l-)(L+l)/L e 2, l p 1 M/L 

is a second degree character associated to £. Applying Theorem II . II now finishes the 
proof. □ 
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As in the discrete case we need to define a special bijective map in order to give 
the definition of a Wilson basis associated with a lattice with arbitrary canonical 
generator matrix. 

Lemma 4.7. Let ^,6 G Z^ ™£/i < b < ^- and let a, (3,m,Q,no : c,d be defined 
as before. Then the mapping ip : Z 2 — > Z 2 defined by 



ip(m, n) 



is bijective and we have 



(mm,Q 



2pd 



(m, n) 



6 = 0, 



{(m mod £,n mod 2c) : (m,n) G ^ _1 ({0, . . . , 2p - 1} x {0, . . . , | - 1})} 



{0,...,£-l}x{0,...,2c-l}. 



Proof. The proof of this lemma is very similar to the proof of Lemma 13.71 we 
therefore omit it. □ 

Note that the following definition reduces to Definition ^. 4l in the case of a diagonal 
canonical generator matrix. 

Definition 4.8. Let A be a lattice in It, x Z^ with canonical generator matrix 
A given by 



A 



2p 

p 



Let g be some L -periodic function on Z, and let ip be defined as in Lemma \4- 7| Then 
the Wilson system W(g, A, C L ) = {V'm,n}(m,n)6/; where L = {0, . . . ,p— 1} x {0, ^} U 
{0, . . . , Ip - 1} x {1, . . . , ^ - 1}, is given by 



V'm.n — 3jj 1 (2m,ii)i+if 2 (2m,n)f),j32(2m,ii)p! if 171 — 0, 

and for m = 0, . . . , 2p — 1, n = 1, . . . , ^ — 1, 

+9 Vl ( m ,-n)±+tp2(m,-n)b,tp2(m,-n)p)' if m + 11 even, 

1 (m,n) ^+iy3 2 (™,n)'),¥'2(m,n)p 
~9ip 1 (m-n)^+i P 2(m-n)b, V 2(m-n)p)i if m + n odd. 



p-l, n = 0, ^, 



The following theorem is the analog to Theorem 13.91 for the space C L . 
Theorem 4.9. Let A be a lattice in Z^ x Z^ with canonical generator matrix A 
given by 



A 



2p 

p 



Let g be some L-periodic function on Z such that U 1 g is real-valued, let M := 2p, 



N 



2c 



, M := 2q, N 



id let U and tp be defined as in Proposition \4-(>] 



id Lemma \4- 7\ respectively. Then the following conditions are equivalent. 
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(i) {s m ^ +n6) „ p }m=0,...,M-i,n=0 1 ...,iV-i is a tight frame for C L with frame bound 
2. 

(ii) {(^ _1 g)„ l ^,„ (? } m= o,...,A7-i,, l= o N-i is a tight frame for C L with frame bound 

2. 

(iii) W(U~ 1 g, x q'El, i C L ) is an orthonormal basis for C L . 

(iv) yV(g,A,C L ) is an orthonormal basis for C L . 

Proof. Let a be defined as before and let tp = ((pi, ip 2 )- Then we compute 
a(cpi(m,n)-^ + ip 2 (m,n)b,(p 2 (m,n)p) 

= (a(ipx(m,ri)^ + ip 2 (m,n)b) + f3(p 2 (m,n)p, 7(951 (m,n)| +ip 2 (m,ri)b) + 5(p 2 (m,ri)p) 
= (m{a±m + abn Q + f3pn ) + n(-a^^ + ab^^ + Pp^^), 
m( 7 Am + 7 6n + Spn ) + n(-7^ + + Sp^)) 

= ( mc > m i ^L h A h m o + bn o) - ^p^o) 



_L_ n ( t d j t bd t pd \\ 

-itTno+bno a Jpmo+bno ab+f3p pno ab+[3p'' 

. „-abpn -f3p 2 no+ubpn -up£ ; m —apbno-.-. 

= (mc,nta( rT — — . , . , „ f )) 

v ' v a(^m +bn )(ub+f3p)pn " 

1 ,11 a£-mo+abno+/3pno Nx 
= (me, —ntdpl . L — — — , , , )) 

y a(^mo+bn )(ab+f3p)pn " 

= (mc. —ntp—-r , C f , 7 ; g . )) 

^ a(^m +bno)(ab+l3p)pn " 



— (mc — nt — ? ~ ) 

' ( 35 m +bn )pn 1 



2p' 



where in the last step we used (|4.5ll and st — —{■§p~ m o + bno)pno. Since \a\ — 1, we 
have c = s = gcd(a^,ab + (3p) is a factor of ^ and hence of -|. Using Lemma l4~71 
the equivalence of (i) and (ii) follows immediately from Proposition 14.61 since U is 
unitary and |C(m, n)\ = 1. Proposition ^. 5l proves (ii) (iii). Therefore it remains to 
prove the equivalence of (iii) and (iv). For this, we will use the following implication 
of Proposition ^. 61 



U(U-'g) m ^ nq = U{U-'g) 

— C(ip(m,n)) 1 5 ipi ( m , n )i+ l/ , 2 ( m ,„)b, ¥)2 ( m ,„)p- 



An easy but tedious calculation shows that C(<p(m, n))^ 1 does not depend on the 
sign of n. Using now the definition of a Wilson basis, the fact that U is a unitary 
operator, and the fact that \C(ip(m, n))| = 1 yields the result. □ 

Tight Gabor frames in C L can be constructed in the same way as for £ 2 (Z) and 
i 2 (R) by using the "inverse square root trick" . Furthermore, it has been shown in |27j 
that for properly localized windows the dual window constructed in C L by "sampling 
and periodization" of the frame {g m a,nb} converges to the dual window S _1 g with 
increasing sampling rate and increasing periodization interval, see |27| for details. 
This result can be easily extended to tight windows. We refer also to [251 1321 ITU] for 
related results and leave the details to the reader. To obtain tight Gabor frames in 
C L that satisfy the required conditions of the theorem above and have good time- 
frequency localization one can thus essentially proceed analogous to the example at 
the end of Section [21 
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5. Conclusion. We have demonstrated that ortho-normal Wilson bases for L 2 (R) 
(with excellent time-frequency localization) can be constructed for general time-fre- 
quency lattices. Of course any numerical implementation has to be done in a discrete 
setting. Somewhat longer proofs establish a similar result for the spaces £ 2 (Z) and 
C L for non-rectangular time-frequency lattices. The approach based on metaplec- 
tic transforms used in this paper suggests that the main results can be extended to 
the setting of symplectic time-frequency lattices on general locally compact abelian 
groups. 

Furthermore, our results imply that from a practical viewpoint it is indeed possible 
to extend OQAM-OFDM or cosine-modulated filter banks to general time-frequency 
lattices. Moreover, we expect that the benefits of using general time-frequency lat- 
tices will be even more pronounced for images and higher-dimensional signals. Our 
expectation is based on the fact that in the theory of sphere packings (and sphere 
coverings) the advantages of the optimal sphere packing over the packing associated 
with the rectangular lattices increases significantly with the dimension of the space [S]. 

An interesting research problem is thus to investigate how to extend the results 
in this paper to L 2 (M. d ) for non-symplectic lattices as well as to find optimal time- 
frequency lattices in R 2d for d > 1. One possibility to define an "optimal" time- 
frequency lattice is to fix the function g to be a Gaussian, say, and then find that 
time-frequency lattice of fixed density which minimizes the condition number of the 
associated Gabor frame operator as indicated in |33|. 

Acknowledgments. We thank the referees for valuable comments and sugges- 
tions which lead to an improvement of the results and presentation in this paper. 
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